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_j- ' Abstract 

In this paper, we investigate the error performance of the protograph (LDPC) codes over Nakagami- 

m fading relay channels. We first calculate the decoding thresholds of the protograph codes over such 

channels with different fading depths (i.e., different values of m) by exploiting the modified protograph 

extrinsic information transfer (PEXIT) algorithm. Furthermore, based on the PEXIT analysis and using 

Gaussian approximation, we derive the bit-error-rate (BER) expressions for the error- free (EF) relaying 
^ ■ 

protocol and decode-and-forward (DF) relaying protocol. We finally compare the threshold with the 

\Q ' theoretical BER and the simulated BER results of the protograph codes. It reveals that the performance of 

\Q '. 

• , DF protocol is approximately the same as that of EF protocol. Moreover, the theoretical BER expressions, 

"xl" . 

which are shown to be reasonably consistent with the decoding thresholds and the simulated BERs, are 



m 



able to evaluate the system performance and predict the decoding threshold with lower complexity as 
compared to the modified PEXIT algorithm. As a result, this work can facilitate the design of the 
S^ ' protograph codes for the wireless communication systems. 

Index Terms 

Extrinsic information transfer (EXIT) algorithm, Gaussian approximation, Nakagami-m fading, Pro- 
tograph low-density parity-check (LDPC) code, Relay channels. 



I. Introduction 

Spatial diversity is an effective technique to enhance the quality and reliability of wireless communi- 
cations and can be typically achieved using multiple antennas at a transmitter and/or a receiver. Besides, 
spatial diversity can also be obtained by the use of relaying. The relay channel, which consists of a source, 
a relay, and a destination, was probably first proposed in the 70s [1]. Later, Cover and Gamal have further 
developed the theory of such a channel [2]. Recently, half-duplex relaying [3] has been verified to be 
relatively more practical in comparison with the full-duplex one due to simpler implementation. 

To improve the performance, forward error correction (FEC) codes have been applied to relay channels. 
As a capacity-approaching code, Low-density parity-check (LDPC) code has been developed to signif- 
icantly enhance the performance of relay channels for different channel conditions, such as non-fading 
[3-5] and fading scenarios [6,7]. Moreover, [8] provided the corresponding analysis of the LDPC-coded 
relay systems. At the same time, there has been a growing interest in protograph LDPC codes. It was 
shown in [9-11] that protograph (LDPC) codes not only achieve superior performance but also possess 
simple structures to realize linear encoding and decoding. For this reason, protograph codes have been 
used in additive white Gaussian noise (AWGN) relay channels [12], partial response channels [13], and 
Rayleigh fading channels [14, 15]. However, thus far, its analytical performance for general fading relay 
channels is not very well understood, which has motivated the results of this paper. 

In this paper, we study the error performance of a half-duplex protograph LDPC-coded relay system 
over Nakagami-m fading channels. The Nakagami-m fading channel is a general type of fading channel 
that encompasses Rayleigh fading (m = 1) and AWGN (m — > oo) channels as special cases. We first use 
the modified protograph extrinsic information transfer (PEXIT) algorithm [15] to analyze the decoding 
threshold of the protograph codes with different fading depths (i.e., different values of m). Afterwards, 
we derive the bit-error-rates (BERs) for the error-free (EF) and the decode-and-forward (DF) relaying 
protocols exploiting the modified PEXIT algorithm and using Gaussian approximation [16], and find that 
the BER analytical method has lower computational complexity than that of the modified PEXIT one. 
Simulated results show that the performance of the DF protocol approaches very close to that of EF, 
which align well with the analytical results. Furthermore, the results also indicate that the performance of 
the codes is improved as the fading depth decreases (higher m), but the rate of improvement is reduced 
simultaneously. 

The remainder of this paper is organized as follows. In Section II, the system model over the Nakagami- 
m fading channel is described. In Section III, we analyze the protograph LDPC coded relay system 



using the modified PEXIT algorithm. Moreover, we derive the BER expressions of the protograph codes. 
Numerical results are carried out in Section IV, and conclusions are given in Section V. 

II. System Model 

We consider a two-hop half-duplex relay system model with one source, S, one relay, R, and one 
destination, D, as shown in Fig. 1. Each transmission period is divided into two time slots, with the first 
slot being the broadcast time slot and the second slot being the cooperative time slot. During the first 
time slot, S broadcasts the codeword to other terminals (including R and D). Then, in the second time 
slot, R cooperates with S to forward the re-encoded message of S to D while S remains idle. 1 During 
each transmission period, D stores the received codeword for decoding at the end of the second time 
slot. Mathematically, the received signals can be written as 

rm,j = hsRjXj + nsn,j, (1) 

rm,j = hsDjXj + n S Dj, (2) 

r-D2,j = h-RD,jXj + n RDj -, (3) 

where, Xj and Xj denote the binary-phase-shift-keying (BPSK) modulated signals corresponding to the 
jth coded bit Vj and re-encoded bit Vj , respectively; trij, rr>i,j, rD2,j represent the received signals 
of the jth coded bit at the relay in the 1st time slot, at the destination in the 1st time slot, and at 
the destination in the 2nd time slot, respectively; ksR,j, hso,j, and ^rdj are the mutually independent 
Nakagami-m fading channel coefficients of the S-R link, the S-D link, and the R-D link, respectively; 
and nsRj, "-sdj, and tirdj are the AWGN with zero mean and variance of a\ {a\ = ^£-). 

We define the channel gain of the S-R link as 7srj = |^sr,j | 2 (similar expressions will be assumed 
for other links). Moreover, we assume that the receiver (R or D) knows perfectly the channel state 
information (CSI) for decoding. At the destination, the signals from the source and relay are combined 
by a maximum ratio combiner (MRC). To simplify the analysis, we also assume that the distance from S 
to D is normalized to unity, while the distances from S to R and from R to D are d and l — d, respectively. 2 
Based on this assumption, we have hsvj = a SD,j, ^SR,j = as ?' 3 , and /irdj = "^ l ' wnere a * s trie 
Nakagami-m fading parameter with the probability density function (PDF) expressed as 

^-^(^"^'-"(-n**)- form£05 ' (4) 

'in our model, we assume that the relay and the source adopt the same coding scheme. 

2 Note that the assumption of collinearity of S, R and D will not affect any of the derivation of our results. 



where m is the fading depth, = E[a 2 ] = 1, E[-] is the expectation operator, and T(-) is the Gamma 
function. In this paper, we focus on ergodic channels, in which the channel fades significantly rapidly 
such that it varies bit by bit. 3 

III. Performance Analysis 

In this section, we derive the BER expression of the protograph code in our system with the EF and 
DF protocols based on the modified PEXIT algorithm and using Gaussian approximation [16]. In the 
analysis, it is assumed that the all-zero codeword is transmitted and the block length of the code is infinite 
[17]. The output (extrinsic and a-posteriori) log-likelihood-ratio (LLR) messages of Nakagami-m fading 
channels are approximately following the symmetric Gaussian distribution [18]. 

Firstly, we review the protograph code [9]. A protograph is a Tanner graph with a relatively small 
number of nodes. A protograph code is a large protograph (called derived graph) obtained by the "copy- 
and-permute" operation on a protograph. Accordingly, a code with different block lengths can be produced 
by different number of times of the "copy-and-permute" operation. A protograph with N variable nodes 
and M check nodes can be described by a base matrix B = (bij) with dimensions M x N, where bij 
denotes the number of edges connecting the variable node Vj to the check node c,. Since the protograph 
code always has some punctured variable nodes, we define Pj as the punctured label of Vj (Pj = if Vj 
is punctured; otherwise Pj = 1). We also define the following LLRs and mutual information (MI). 

LLR value: 

• Lxv{i,j) denotes the input (a-priori) LLR value of Vj corresponding to Cj on each of the bij edges. 

• Lxc(i,j) denotes the input (a-priori) LLR value of Cj corresponding to Vj on each of the fry edges. 

• Lft v (i,j) denotes the output (extrinsic) LLR value passing from Vj to the c%. 

• L/Ecihj) denotes the output (extrinsic) LLR value passing from a to the Vj. 

• ^app(j) denotes the a-posteriori LLR value of Vj. 
MI: 

• I& y (i,j) denotes the a-priori MI between Lj± v (i,j) and the corresponding coded bit Vj. 

• I^dhj) denotes the a-priori MI between L\ c (i,j) and the corresponding coded bit Vj. 

• lEv{i,j) denotes the extrinsic MI between Le v (z, j) and the corresponding coded bit Vj. 

• lEc(i,j) denotes the extrinsic MI between Le c (^,j) and the corresponding coded bit Vj. 

3 Although we consider the ergodic (i.e., fast fading) scenario here, the results in this paper are also applicable for the 
quasi-static case (where the fading parameter of each link is kept to be constant for a code block). 



• ^app(j) denotes the a-posteriori MI between L app (j) and the corresponding coded bit Vj. 

Note that the subscripts "SR", "SD", "RD", and "D" are used to denote the S-R link, the S-D link, 
the R-D link, and the destination, respectively. For example, Lo,R-Av(hJ) represents the a-priori LLR 
value of Vj corresponding to a on each of the b%j edges of the S-R link. 

A. BER of EF Relaying Protocol 

At the destination receiver, the output of MRC corresponding to the jth coded bit is given by 

ynj = hsDjroij + hnr>,jrm,j- (5) 

Afterwards, the initial LLR value ^D-ch,? of the jth coded bit is calculated as 

Pi(v j =0\y Dd ,h j ) Pr(xj = +l|ltoj, hj) 

-^D-chj- = m ^—< : — r = in ■ 



Pt(vj = l\y D j,hj) Pt(xj = -1\vd,j, hj) 

*VDj _. 2 



<?l cl 



(hsDjrmj + hKDjr T) 2,j} 



— [hsn,j(hsT>,jXj + nsD,j) + /ird,j(/ird,j£j + ^rdj)] 

— [\hsn,j\ Xj + |/irdj| Xj + hsD,jnsT>,j + hBD,jnRD,j] 



0", 

2 



[(TSDj + 7KD,j)xj + hsDjnsDj + /iRDJ^RDj] , (6) 



where Pr(-) denotes the probability function, h_, = {hsR,j,hsD,j, ^RDj} is the channel vector, 7sdj = 
I^SDjI 2 = "Id ?' 7RD,j = I^RDjI 2 = (i^d) 2 ' anc * x i = % ( err or-free decoding at the relay). 

Based on the all-zero codeword assumption (xj = +1), given a fixed channel realization (a fixed fading 
vector hj), we can easily obtain the expectation and the variance of Ld-cIij, resulting in 

2 2 

IE[L D _ chiJ -] = — (tsdj + 7rd,j) = — Ad j, (7) 

var[L D -chj] = — (7SDjc n + 7RD,jO" n ) = — Adj, (8) 

where Ad j is the short-hand notation of 7sdj + 7rd j- As seen from (7) and (8), Ld-cIiJ follows the 
symmetric Gaussian distribution, i.e., -Lo-chj ~ A^-^Adj, -4-Adj). Subsequently, substituting o~\ = 
2R(E' In ) ^ nto $) (-^6 = 2 ' tne normalized factor 1/2 is used to keep the total energy per transmitted 
symbol to be E s ) and considering the punctured label Pa, we have 

var[L D _ chj ] = SRPjXvjiE'jNo) = ARPjX^Ei/No). (9) 



Thus, the modified PEXIT algorithm [15] can be applied to our system by using the expression (9). 4 
For the qth (q = 1, 2 . . . , Q, where Q is the total number of channel realizations 5 ) channel realization 



h q j, the a-posteriori LLR of Vj during the (t + l)th iteration is written as 

M 
^D + -app >(? (i) = E M^Av,g&i) + L D-ch, g J. (10) 

The expected value of L^L Av q {hj) is then expressed by 

1 Q 

l£\ pp (j) = nL^ m (j)} = ^E L D + - apPl9 (i)- (id 

As the output LLR values, i.e., the extrinsic LLR and a-posteriori LLR, are approximated to follow a 
symmetric Gaussian distribution, we can evaluate the variance of L^_ (j) as 

where /^(j) = l^app,^)] = £ £?=i #4^0') « *e expected value of lg^V) over all 
channel realizations, the expression of the a-posteriori MI i^L aP p „(j) is obtained as [15, (21)], denoted 
b y J D + - ap p, ? 0') = HId\ v ^J), ■ ■ ■,lg~- Av , q (M,j);vax[L B „ ch>q;j }), and J(a ch ) is the MI between a 
BPSK modulated bit and its LLR value L cn ~ J\f (^f- , a^J over an AWGN channel, represented as [17] 

J(a ch ) = l- V 7 log 2 [1 + exp(-e)] de (13) 

The corresponding inverse function is given by [17] 

_, 771X 2 + r]2X + rjz^fx if < x < 0.3646, 

J x (x) = < (14) 

774 In [775(1 — x)l + r/6X otherwise, 

where 771 = 1.09542, 772 = 0.214217, 773 = 2.33737, 774 = -0.706692, 775 = 0.386013 and t? 6 = 1.75017. 

With the help of (12)— (14), the BER of the jth variable node after t iterations is evaluated by 

P£-\(j) = \ Pr^-appti) < 0|x, = +1) + \ Pr(L*+! app (i) > 0|xj = -1) 




E^iappO - )!^ = +1] \ 1 , / var^+Wi)!^ = +l]/2 



Zvar^^pOOl^ = +1] / l \ j2var[Z^ 1 app (j)|x J = +1] 



var[^ + _ap P (j)l^ = +l]\ 1 e {J-Hrg'ti)) 



2\/2 I 2 V 2^2 




(15) 



4 According to [15], the maximum iteration number of the modified PEXIT algorithm T^ax should be large enough in order 
to ensure the complete convergence of the decoder, i.e., T^ax > 500. 

5 To ensure the accuracy of our derived BER expressions, Q should be set to a sufficiently large integer, i.e., Q > 10 5 . 



where erfc(-) is the complementary error function, defined as erfc(x) = -j= J^° e T "dr. 
Finally, the averaged BER of a protograph code after t iterations is written as 



1 ,rrc °" 1 ^^' 
3=1 "' i=i 



N N 

pt+i _ pt+l - _!_ V P t+1 ( i) - — V erfc 



iV^ D ~ 6W 2iV^ I 2^2 



(16) 



£. BER of DF Relaying Protocol 

In DF, if the relay can decode received signals correctly, it re-encodes the decoded information and 
then forwards them to the destination; otherwise, it does not send message or remains idle. Consequently, 
the BER of the DF protocol of the jth variable node after t iterations can be found as 

p&Uj) = P&- b (J)P&UJ) + [i - r&M Pht\U), (iv) 

where P^_ b (j), ^gD-ftC?)' anc ^ P^-biJ) are trie corresponding BERs of the jth variable node after t 
iterations at the relay receiver with the signal from source, at the destination receiver with one signal 
from the source, and at the destination receiver with two signals both from the source and relay (i.e., 
BER of the EF protocol), respectively. 

The initial LLR of the jth variable node of the S-R link is expressed by [18] 
T 2/isRjmij 2/isRj, 2 

In (18), we have Ef-LsR-chj] = 7 ^?' J and var[LsR-ch,j] = 7 ^?' J ■ Hence, during the (t+l)th iteration, the 
a-posterior MI of the jth variable for the gth (q = 1,2, . . . ,Q) channel realization and the corresponding 
expected value are respectively given by 

'sR-apP,^') = P^R-A^J), ■ ■ • ,4R-Av )9 (^i);var[L S R-ch, gj ]), (19) 

1 Q 
^app(i) = E[^app, f/ (i)] = qE ^-app^O')" (20) 

g=l 

Further, the variance of the expected a-posterior LLR associated with (20) is yielded in terms of the 
Gaussian assumption 

var^sJtappW] = {^SR-appO-))} 2 ■ (2D 

Note that E[L^ app (j)] = var[L i + 1 _ app (j)]/2. 

Using (21), we can therefore obtain the BER of the jth variable after t iterations as 



><+l f^ _ i erfc ^SR-appU^J ~ t*] = l erfc o I'SR-ap p 

2 I . /wrS mi*.- =-1-11/ 2 I 2^2 



PsZUti) = o erfc i oa -™~" = | = ^erfc [ - "Tr ] • (22) 



Likewise, the BER of the jth variable after t iterations of the S-D link can be written as 



yt+l (a\ - l or f„ I ^SD-a pp 



p&uti) = o erfc — : ,r I (23) 



with the parameters subjected to 

1 Q 

-^SD-appO) = E [^SD-app,qU)] = 7> / ^SD-app.ou) 



(24) 



r r 1 4 7SD,g J 

where LsD-ch,gj is the initial LLR of Vj for the qth S-D link realization and 7SD,g,j = |^SD,g,j| 2 - 

Substituting (15), (22) and (23) into (17), the BER of the jth variable node with DF after t iterations 
can be formulated. Finally, we get the averaged BER with DF protocol after t iterations as 

1 N 

3=1 

Note also that 

• The maximum number of iterations of the theoretical BER analysis T max (t = T max ) equals to that 
of the simulations (T max is always much smaller than T^ ax ). In this paper, we set T max = 100 as 
in [15,20]. Therefore, its computational complexity can be reduced by v ma ^, P max; as compared to 
the modified PEXIT algorithm. 

• The BER analysis can be used to evaluate the performance of the protograph codes for any E^/Nq 
and T max , while the modified PEXIT algorithm can only be exploited to derive the E^/Nq threshold 
above which an arbitrarily small BER can be achieved (i.e., all the code blocks can be successfully 
decoded) for a sufficiently large T^ 



■max' 



IV. Numerical Results and Discussions 

In this section, we firstly analyze the decoding threshold of two typical protograph codes, i.e., the AR3A 
and AR4JA codes, over Nakagami-m fading relay channels utilizing the modified PEXIT algorithm [15]. 
Then, we compare the decoding thresholds, the theoretical BER and simulated BER results of the two 
protograph codes. For all the following results, the distance of S-R link d is set to 0.4. 

The AR3A and AR4JA codes, which have been proposed by Jet Propulsion Laboratory [10, 11], can 
respectively accomplish excellent performance in the low SNR region and the high SNR region over the 



AWGN channel. The corresponding base matrices of these two codes with a code rate of R = ^±1, 
denoted by Ba3 and Ba4, respectively, are given by 

/ 



2/i 



BA3 



121000 ••• 

021112 1 ••• 2 1 

012111 2 ••• 1 2 

In 



(26) 



B 



A\ 



120000 

031113 1 

yo 1 22 1 1 3 




3 1 
13/ 



(27) 



In (26) and (27), the jth column corresponds to the jth variable node and the ith row refers to the zth 
check node. The variable nodes corresponding to the second columns of the two matrices are punctured. 

A. Decoding Threshold Analysis 

We calculate the decoding thresholds of the AR3A and AR4JA codes with different code rates and 
different fading depths using the modified PEXIT algorithm [15] and show the result in Table I. 6 As 
seen from this table, the threshold of the AR3A code is lower than that of the AR4JA code for a fixed 
code rate and fading depth (fixed R and m). For instance, the thresholds of the AR3A code and the 
AR4JA code are 0.575 dB and 0.722 dB, respectively, with parameters R = 4/5 and m = 2. Moreover, 
the threshold is reduced as the fading depth decreases (higher m) for both the two codes and hence the 
error performance should be improved in the waterfall region. However, the decrease of the threshold is 
reduced when m becomes larger, indicating that the rate of performance improvement is reduced. 

B. BER Performance 

In the sequel, we present the simulated results of the AR3A and AR4JA codes with a code rate of 
R = 4/5 and compare their simulated and theoretical BER curves with the decoding thresholds. We 
denote N p , K p , and L p as the block length, the information length, and the punctured length of the code, 
respectively. Unless specified otherwise, in the simulations, it has been assumed that the decoder performs 
a maximum of 100 iterations for each code block with parameters [N p ,K p ,L p ] = [5632,4096,512]. 

In Fig. 2, we show the BER results (vs. E^/Nq) of the AR4JA and AR3A codes over Nakagami-fading 
relay channels with DF and EF protocols. The fading depth m is set to 2. As seen, the error performance 



6 For the threshold analysis, we only consider the perfect S-R channel (i.e., EF), as is typical in many practical cases [19]. 
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of the codes with the DF protocol approaches closely to that with EF, which suggests that the relay 
can decode most of the received code blocks successfully. Moreover, the AR3A code outperforms the 
AR4JA codes for the range of E^/Nq under study. At a BER of 10~ 5 , the AR3A code achieves a gain 
approximately 0.2 dB as compared to the AR4JA code both for the DF and EF protocols. In the same 
figure, consider the AR3A code with DF protocol at a BER of 10~ 5 , the theoretical BER result is in 
good agreement with the PEXIT threshold within 0.2 dB. Moreover, the corresponding simulated curve 
has another gap about 0.7 dB to the theoretical one, which is reasonably consistent with the results in 
[16,20]. It is because our theoretical BER formulas are derived subjecting to the infinite-block length 
assumption. 

Fig. 3 presents the BER results (vs. E^/Nq) of the AR3A code in the EF relay systems for different 
fading depths (m = 1, 2, 3, and 4). As E^/Nq increases steadily, both the two codes start to perform better 
for a larger value of m in terms of the thresholds, the theoretical and simulated BER curves. Nevertheless, 
the improved gain is reduced as m increases. For example, at a BER of 10 -5 , the improved gains are 
about 1.2 dB, 0.45 dB, and 0.25 dB as m increases from 1 to 2, 2 to 3, and 3 to 4, respectively. Moreover, 
it can be observed that gaps between the theoretical BERs and the simulated results are around 0.8 dB, 
0.7 dB, 0.55 dB, and 0.5 dB for m = 1, 2, 3, and 4, respectively. Simulations have also been performed 
for the AR4JA code and with the DF protocol, and similar observations are obtained. 

V. Conclusions 

The performance of the protograph LDPC codes over Nakagami-m fading relay channels has been 
studied. The BER expressions for the protograph codes with DF and EF relaying protocols have been 
derived using the modified PEXIT algorithm and Gaussian approximation. The decoding threshold, 
the theoretical BER and the simulated BER results have shown that the error performance of the DF 
protocol is very close to that of EF, which suggests that the relay can decode most received codewords 
correctly. The differences between the threshold and the theoretical BER and between the theoretical 
BER and the simulated BER have been found to be around 0.15 — 0.2 dB and 0.5 — 0.8 dB, respectively, 
showing a reasonable consistence. Consequently, our analytical expressions not only can provide a good 
approximation of the system performance for a large-block length but also predict accurately the decoding 
threshold more efficiently in comparison with the modified PEXIT algorithm. 
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Fig. 1, The system model of the protograph LDPC-coded relay system over fading channels. 
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Fig. 2. BER results of AR3A and AR4JA codes over Nakagami-fading relay channels with the fading depth m = 2. 




Fig. 3. BER results of the AR3A code over Nakagami-fading relay channels with different fading depths (m = 1, 2, 3, and 4). 
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TABLE I 



Decoding thresholds (E b /N ) th (dB) of the AR3A code and AR4JA code with different code rates over 
Nakagami-fading relay channels with different fading depths m = 1, 2, 3, AND 4. 



Code Rate 


AR3A code 


AR4JA code 


m — 1 


m = 2 


771 = 3 


m = 4 


m = 1 


m = 2 


m = 3 


m — 4 


1/2 (n = 0) 


-1.345 


-1.825 


-1.983 


-2.056 


-1.187 


-1.675 


-1.828 


-1.895 


2/3 {n = 1) 


0.004 


-0.759 


-0.994 


-1.098 


0.162 


-0.566 


-0.823 


-0.918 


3/4 {n = 2) 


0.890 


-0.016 


-0.316 


-0.472 


1.128 


0.185 


-0.136 


-0.292 


4/5 {n = 3) 


1.639 


0.575 


0.216 


0.037 


1.805 


0.722 


0.382 


0.192 


5/6 {n = 4) 


2.235 


1.042 


0.653 


0.454 


2.392 


1.164 


0.785 


0.576 


6/7 (n = 5) 


2.732 


1.404 


0.991 


0.778 


2.855 


1.533 


1.092 


0.887 


7/8 (n = 6) 


3.158 


1.748 


1.278 


1.044 


3.284 


1.852 


1.379 


1.150 



